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■ Abstract: We examine the prospects of detecting the hght Higgs boson h'^ of the Constrained 

'"q^. MSSM at the Tevatron. To this end we explore the CMSSM parameter space with > 0, using a 

' , Markov Chain Monte Carlo technique, and apply all relevant collider and cosmological constraints 

D , including their uncertainties, as well as those of the Standard Model parameters. Taking 50 GeV < 

' "ii/2,wo < 4 TeV, \Aq\ < 7 TeV and 2 < tan/3 < 62 as flat priors and using the formahsm of 

Bayesian statistics we find that the 68% posterior probability region for the /i° mass lies between 

X ■ 115.4 GeV and 120.4 GeV. Otherwise, /i" is very similar to the Standard Model Higgs boson. 

5— ( ' 

' Nevertheless, we point out some enhancements in its couplings to bottom and tau pairs, ranging 

from a few per cent in most of the CMSSM parameter space, up to several per cent in the favored 
region of tan P ^ 50 and the pseudoscalar Higgs mass of to^ < 1 TeV. We also find that the other 
Higgs bosons are typically heavier, although not necessarily much heavier. For values of h*^ mass 
within the 95% probability range as determined by our analysis, a 95% CL exclusion limit can 
be set with about 2fb~^ of integrated luminosity per experiment, or else with 4fb^^ (12 fb^^) a 
3cr evidence {5a discovery) will be guaranteed. We also emphasize that the alternative statistical 
measure of the mean quality-of-fit favors a somewhat lower Higgs mass range; this implies even 
more optimistic prospects for the CMSSM light Higgs search than the more conservative Bayesian 
approach. In conclusion, for the above CMSSM parameter ranges, especially mp, either some 
evidence will be found at the Tevatron for the light Higgs boson or, at a high confidence level, the 
CMSSM will be ruled out. 
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1. Introduction 

One of the main goals of the Tevatron and the LHC experimental programmes is to detect 
a Higgs boson. In contrast to the Standard Model (SM), in models with softly broken 
low energy supersymmetry (SUSY), the mass of one Higgs boson is restricted to be fairly 
low, rrih < 150 GeV,^ which allows for a more focused search. On the other hand, even 
in the Minimal Supersymmetric Standard Model (MSSM) there are as many as five Higgs 
mass states: two scalars, and H^, a pseudoscalar, and a pair of charged bosons, 
H^, which makes the experimental search more involved. While Higgs boson tree-level 
mass parameters obey some well-known relations, top and stop loop-dominated radiative 
corrections introduce large modifications to Higgs masses and couplings in terms of several 
unknown SUSY parameters. In the general MSSM most soft SUSY-breaking parameters 
remain fairly unrestricted, which makes it difficult to conduct a thorough exploration of the 
parameter space. Instead, many studies in the general MSSM, including the Higgs discovery 
potential at the Tevatron, often adopt rather arbitrary choices for MSSM parameter values. 

It is therefore interesting and worthwhile to assess Higgs observability in more con- 
strained and well-motivated low energy supersymmetric models. One particularly popular 
framework is the Constrained MSSM (CMSSM), introduced in ref. [^,^ which is defined 
in terms of the usual four free parameters: the ratio of Higgs vacuum expectation values 
tan /3, as well as the common soft SUSY-breaking parameters for gauginos mi/2, scalars mg 
and tri-linear couplings ^o- The parameters w-o and are specified at the GUT 

scale Mgut ^ 2 x lO^^ GeV and serve as boundary conditions for evolving the MSSM 
Renormalization Group Equations (RGEs) down to a low energy scale MgusY = ^y^Ji 

^For recent extensive reviews and further references see, e.g., ||^, ^. 

^One well-known implementation of the CMSSM is the minimal supergravity model M. 
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(where mj-^ denote the masses of the scalar partners of the top quark), chosen so as to 
minimize higher order loop corrections. At MgusY the conditions of electroweak symmetry 
breaking (EWSB) are imposed and the SUSY spectrum is computed at mz- The sign of 
the Higgs/higgsino mass parameter /i, however, remains undetermined. 

Prospects for Higgs collider searches in the CMSSM and other unified models have 
been explored in several recent analyses [||, A usual approach is to perform a fixed 
grid ( "frequentist" ) scan in some of the CMSSM parameters (typically and mo) while 
keeping the remaining ones (typically tan/3 and Aq) and also SM parameters fixed. The 
resulting "allowed regions" of parameter space then often largely underestimate the true 
extent of the uncertainties, mainly because of the existence of degeneracies in parameter 
space that this procedure does not account for. On the other hand, a full scan over a 
parameter space of even moderate dimensionality using grid techniques is highly inefficient. 
Not only the size of (and time spent on) the scan grows as a power-law with each new 
parameter added to the scan, but there are other limitations. It is difficult to incorporate 
residual error-bars of relevant SM parameters, which are often simply fixed at their central 
values. Experimental limits on SUSY are applied at some arbitrary confidence level, e.g., 
at 1 or 2(7. As a result, it is very difficult, if not impossible, to derive a global picture of 
the most probable ranges of SUSY parameters. 

Recently more efficient exploration methods based on the Markov Chain Monte Carlo 
(MCMC) technique have been successfully applied to studying SUSY phenomenology 
and are becoming increasing ly popular [|, |, 0, |ll|, |l| . The MCMC technique allows one to 
make a thorough scan of a model's full multi-dimensional parameter space. Additionally, 
by combining the MCMC algorithm with the formalism of Bayesian statistics, maps of 
probability distributions can be drawn not only for the model's parameters but also for all 
the observables (and their combinations) included in the analysis. 

In our first paper we applied this approach to performing a full analysis of the 
CMSSM. As in a similar (and concurrent, but independent) work of Allanach and Lester 
we applied all relevant constraints on Higgs and superpartner masses from collider searches, 
from the rare processes BR{B — > ^^7) and BR{Bs — > from the anomalous mag- 

netic moment of the muon (g — 2)^, and also from cosmology, on the relic abundance of the 
lightest neutralino O^/i^ assumed, in the presence of i?-parity, to be the cold dark matter 
in the Universe. We also took into account residual error bars in the pole top mass, the 
bottom mass and as{Mz)^^^ ■ Going beyond the work of ref. |p, we further included in our 
analysis the experimental error in the fine structure constant measured at mz (which had 
a sizable impact on Jl^/i^), explored wider ranges of mo up to 4 TeV (which allowed us to 
explore the focus point region) and computed the spin-independent cross section for dark 
matter neutralino scattering off nuclei dp ^ (but did not use it as a constraint on the model 
because of astrophysical uncertainties) . We further included constraints from contributions 
to myy and sin^ O^s from full one-loop SM and MSSM corrections and from two-loop SM 
corrections involving the top Yukawa. Last but not least, we also emphasized the difference 
between high posterior probability regions of the parameters in the Bayesian language and 
those of the high mean quality-of-fit, i.e., possibly limited ranges of parameters that give 
the best fit to the data. We have found that in the CMSSM the two can be rather differ- 
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ent, which is a consequence of the complex dependence of the model's parameters on the 
applied constraints. The main results of both analyses and |11] came out remarkably 
consistent with each other, in spite of the above differences and, additionally, of some nu- 
ances in computing the likelihood function. In particular, high probability regions showed 
preference for mi/2,n^o ^ 1 TeV, but not for values nearly as low as claimed in ref. |13| 
based on a analysis. 

In our present work we apply an important recent shift in the SM value of BR{B 
Xs7). In 1 11] we used the previous SM prediction (3.70 ±0.30) x 10""^ HI, which included 
a full NLO calculation and partial charm mass contribution. Recently, partial NNLO 
contributions, most importantly an approximate charm mass one, have been obtained 
in [^] which led to a rather dramatic shift down to (3.15 it 0.26) x 10~^, with a further 
slight decrease to (2.98 it 0.26) x 10~^ after including some additional subtle effects due to a 
treatment of a photon energy cut > 1.6 GeV [jl6|. At the same time, the experimental 
world average has somewhat increased from (3.391q27) ^ 10~^ to the current range of 
(3.55 ± 0.26) X 10~^ ||l^. This leads to some discrepancy, at the level of lAa, between 
the SM and the experimental average. More importantly, because the SM value for the 
branching ratio has moved down below the measured one, any potential overall SUSY 
contributions should now preferably be constructive, in contrast to the situation before. 
This will lead to significant shift in the preferred regions of the CMSSM parameter space. 

In the present work we further include full two-loop and available higher order SM 
corrections, as well as dominant two-loop MSSM gluon corrections, to mw and sin^ 0eff- 
Assuming for comparison the previous values of BR{B Xsj), including these observ- 
ables does not, however, lead to any appreciable differences with respect to [||, |l^], as was 
also shown in very recent updated analysis of Allanach et al., |jl^ who included them at 
comparable level. We also update several experimental constraints, as discussed below. 
We additionally compute Bg mixing, AM^^, which has recently been precisely measured 
at the Tevatron by the CDF Collaboration |18]. 

We devote this paper to a study of the light Higgs boson in the CMSSM and to the 
prospects for its detection at the Tevatron. Results of our new analysis regarding the 
CMSSM parameters will be presented elsewhere |19]. 

Our main results can be summarized as follows. We find that, imposing flat priors 
on wide ranges of the CMSSM parameters: 50 GeV < mi/2,mo < 4 TeV, \Ao\ < 7 TeV 
and 2 < tan /? < 62, the 68% probability region for the mass of the lightest Higgs is given 
by 115.4 GeV < ruh < 120.4 GeV, the 95% probability range being 112.5 GeV < rrih < 
121.9 GeV. Its couplings generally closely match those of the SM Higgs boson with the 
same mass, although we find some differences at the level of a few to several per cent. 
Ensuing prospects for experimental Higgs search at the Tevatron look excellent. So far, 
with about 1 fb~^ of data analyzed, both CDF and DO Collaborations have been able to 
put interesting limits on the Higgs cross sections [20| for some specific choices of the general 
MSSM parameters which are, however, not representative of unified models. As 

more data are coming in, both collaborations will soon be in a position to start probing 
unification-based models, including the CMSSM. In the whole Higgs mass range given 
above, a 95% CL exclusion limit on the SM-like Higgs boson can be set with about 2fb~^ 
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of integrated luminosity per experiment. In the CMSSM a 3a (5cr) signal should be seen 
in this mass range with about 4fb"^ (12 fb"^) of data per experiment. On the other 
hand, with about 8fb^^ of integrated luminosity eventually expected per experiment, a 5a 
discovery wih be possible, should the light Higgs mass be around 115 GeV. In conclusion, if 
the CMSSM (or another supersymmetric model with similar light Higgs boson properties) 
has been chosen by nature, then the Higgs boson with SM-like properties will be discovered 
at the Tevatron. 

The paper is organized as follows. In section |2| we briefly summarize the main features 
of our Bayesian analysis and then provide our updated list of experimental constraints. 
We then proceed to present, in section ^, our results for Higgs mass distribution and other 
properties. Section ^ is devoted to a discussion of light Higgs production and decay at the 
Tevatron. In section ^ we present our summary and conclusions. 



2. An Outline of the Phenomenological Analysis 
2.1 Posterior probabilities 

Our procedure based on MCMC scans and Bayesian analysis has been presented in detail 
in |ll|. Here, for completeness, we summarize its main features. 

We are interested in delineating high probability regions of the CMSSM parameters. 
We fix sign(/i) = +1 throughout and denote the remaining four free CMSSM parameters 
by the set 

6* = (mo,mi/2, Ao,tan/3). (2.1) 



As demonstrated in |^, 11], the values of the relevant SM parameters can strongly influence 
some of the CMSSM predictions, and, in contrast to common practice, should not be 
just kept fixed at their central values. We thus introduce a set ^jJ of so-called "nuisance 
parameters". Those most relevant to our analysis are 

= iMt,mbimbf^,aemiMzf^,as{Mzf^), (2.2) 

where Mt is the pole top quark mass. The other three parameters: mh{mb)^^ — the 
bottom quark mass at mi,, aem(-^z)^^'^ and as{Mz)'^^ — respectively the electromagnetic 
and the strong coupling constants at the Z pole mass Mz - are all computed in the MS 
scheme. 

The set of parameters 6 and ip form an 8-dimensional set m of our "basis parameters" 
m = {9,'tp).'^ In terms of the basis parameters we compute a number of collider and 
cosmological observables, which we call "derived variables" and which we collectively denote 
by the set ^ = (^i, ^2, • • •)• The observables, which are listed below, will be used to compare 
CMSSM predictions with a set of experimental data d, which is currently available either 
in the form of positive measurements or as limits. 

In order to map out high probability regions of the CMSSM, we compute the posterior 
probability density functions (pdf's) p{m\d) for the basis parameters m and for several 



^In we denoted our basis parameters with a symbol 77. 
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SM (nuisance) parameter 


Mean value Uncertainty 
^ a (exper.) 


Ref. 


Mt 

l/aem(Mz)^^ 


171.4 GeV 2.1 GeV 
4.20 GeV 0.07 GeV 
0.1176 0.002 
127.955 0.018 


P 

[24 
[24 
[24 





Table 1: Experimental mean /i and standard deviation a adopted for the likelihood function for 
SM (nuisance) parameters, assumed to be described by a Gaussian distribution. 



observables. The posterior pdf represents our state of knowledge about the parameters m 
after we have taken the data into consideration (hence the name). Using Bayes' theorem, 
the posterior pdf is given by 

On the r.h.s. of eq. ( |2.3D , the quantity p{d\£^), taken as a function of d for a given m, 
and hence a given ^(m), is called a "sampling distribution". It represents the probability 
of reproducing the data d for a fixed value of £,{m). Considered instead as a function 
of ^ for fixed data d, p{d\(^) is called the likelihood (where the dependence of ^ on m is 
understood). The likelihood supplies the information provided by the data and, for the 



purpose of our analysis, it is constructed in Sec. 3.1 of ref. [11[. The quantity 7r(m) denotes 
a prior probability density function (hereafter called simply a prior) which encodes our 
state of knowledge about the values of the parameters in m before we see the data. The 
state of knowledge is then updated to the posterior via the likelihood. Finally, the quantity 
in the denominator is called evidence or model likelihood. In the context of this analysis it 
is only a normalization constant, independent of m, and therefore will be dropped in the 
following. 

As in ref. [0], 

our posterior pdf 's presented below will be normalized to their maximum 
values, and not in such a way as to give a total probability of 1. Accordingly we will use 
the name of a "relative posterior pdf" , or simply of "relative probability density" . 

2.2 Likelihood function and constraints 



We scan over very wide ranges of CMSSM parameters; compare table 1 of ref. [11[. In 
particular we take flat priors on the ranges 50 GeV < mi/2,mo < 4 TeV, \Aq\ < 7 TeV 



and 2 < tan/3 < 62. (In ref. [11[ we called this the "4 TeV range".) For the SM (nuisance) 



parameters, we adopt a Gaussian likelihood with mean and standard deviation as given in 
table |l|, and we assume flat priors over wide ranges of their values [11[. Note that, with 



respect to ref. [11[, we have updated the values of all the parameters, including the recent 



shift in Mf based on Tevatron's Run-II 1 fb^^ of data. 

The experimental values of the collider and cosmological observables (our derived vari- 
ables) are listed in table |2| and in table 4 of ref. [0], with updates where applicable. In 
particular, in addition to (most importantly) BR{B Xs'y) summarized above, we update 
an experimental constraint from the anomalous magnetic moment of the muon [g — 2)^ 



- 5 - 



Observable 


Mean value Uncertainties 

/i a (exper.) r (theor.) 


Ref. 


71 /f 

sm OeS 

BR(B X,-f) X 10^ 

O ^2 

X 


oU.o9z GeV 29 MeV 15 MeV 
(J.2oioo io X iU io X iU 
28 8.1 1 
3.55 0.26 0.21 
17.33 0.12 4.8 
0.104 0.009 {).in^h'^ 


iTTEIi 

|25] 
|24] 
117] 
|18] 

\f3 




Limit (95% CL) r (theor.) 


Ref. 


sparticle masses 


< 1.0 X 10"^ 

> 114.4 GeV (91.0 GeV) 

f{mh) 

See table 4 in ref. |11 


14% 
3 GeV 
negligible 

I- 


|28] 
|28] 



Table 2: Summary of the observables used in the analysis. Upper part: Observables for which a 
positive measurement has been made. i5a^^^^ denotes the discrepancy between the SM prediction 
and the experimental value of the anomalous magnetic moment of the muon {g — 2)^. For central 
values of the SM input parameters used here, the SM value of BR{B — > Xg^) is 3.11 x 10^^, while 
the theoretical error of 0.21 x 10"'' includes uncertainties other than the parametric dependence 
on the SM nuisance parameters, especially Mt and as{Mz)^^^ ■ As explained in the text, for each 
quantity we use a likelihood function with mean ^ and standard deviation s = + t^, where 
a is the experimental uncertainty and r represents our estimate of the theoretical uncertainty. 
Lower part: Observables for which only limits currently exist. The likelihood function has been 
constructed as in ref. ||ll[| , including in particular a smearing out of experimental errors and limits 
to include an appropriate theoretical uncertainty in the observables. 



(denoted here by (5o^^^^) for which a discrepancy between measurement and SM predic- 
tions (based on e+e" data) persists at the level of 2-3cr |24]. We note here, however, 
that the impact of this (somewhat controversial) constraint on our findings will be rather 



limited. We also apply the new values for the measured branching ratio for 6 — > S7 |17| 



and an improved 95% CL limit BR{Bs —>■ fi ) < 1.0 x 10 [27|. In constraining the 



relic abundance ^x^"^ °^ lightest neutralino we use the 3-year data from WMAP |p6|| . 
As a new constraint, we add a recent value of Bg-Bg mixing, AMb^, which has recently 



been precisely measured at the Tevatron by the CDF Collaboration |18]. In both cases 
we use expressions from ref. ||2^ which include dominant large tan /^-enhanced beyond-LO 
SUSY contributions from Higgs penguin diagrams. Unfortunately, theoretical uncertain- 
ties, especially in lattice evaluations of /s^ are still very large (as reflected in table ^ in the 
estimated theoretical error for AMb^), which makes the impact of this precise measurement 
on constraining SUSY parameter space somewhat limited.^ 

For all the quantities for which positive measurements have been made (as listed in 
the upper part of table ^ , we assume a Gaussian likelihood function with a variance given 
by the sum of the theoretical and experimental variances, as motivated by eq. (3.3) in 



*0n the other hand, in the MSSM with general flavor mixing the bound horn AMb^ is in many cases 
much more constraining than from other rare processes pCl. 
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ref. 1 11]. For the observables for which only lower or upper hmits are available (as listed in 
the bottom part of table ^) we use a smoothed-out version of the likelihood function that 
accounts for the theoretical error in the computation of the observable, see eq. (3.5) and 
fig. 1 in lll|. 

The likelihood function for the CMSSM light Higgs boson requires a more refined 
treatment. The final LEP-II lower bound of 114.4 GeV (95% CL) [^] is applicable for the 
case of the SM Higgs boson. It also applies to the lightest Higgs boson h of the MSSM 
when its coupling to the Z boson is SM-like, i.e., when (^^ = g'^ {hZ Z)mssm/ g"^ [hZ Z)sm — 1 
which holds in a so-called decoupling regime of tjia ^ ttt-z- For arbitrary values of niA, 
the LEP-II Collaboration has set 95% CL bounds on ruh and uia as a function of p8| , 
with the lower bound of rrih > 91 GeV for rrifi ~ and <C 1 [^^. In this case we use 
a cubic spline to interpolate between selected points in ruh and translate the above bound 
into the corresponding 95% CL bound in the {mh^Ch) pla^is. We then add a theoretical 
uncertainty Tirrih) = 3 GeV, following eq. (3.5) in ref. |]ll|. (Notice that the parametric 
uncertainties coming from the errors in top quark mass and the strong coupling constant 
have already been fully accounted for by including them as nuisance parameters.) We then 
simultaneously constrain the values of {rrihXh) obtained in the CMSSM by comparing 
them with the 2-dim likelihood function for these two variables from the LEP results. 
Since we find ~ 1 with very high accuracy basically everywhere in the parameter space 
(as we will see later), introducing an extra theoretical uncertainty in (which could be 
implemented by extending eq. (3.5) in ref. [0] to a 2-dim case) would not affect our results 
in any appreciable way.^ This procedure results in a conservative likelihood function for 
nifi, which does not simply cut away points below the 95% CL limit of LEP-II, but instead 
assigns to them a lower probability that gradually goes to zero for lower masses. 

As mentioned in the Introduction, here we make some additional improvements in 
our treatment of the radiative corrections to the electroweak observables Mw and sin^ ^gff- 
We now include full two-loop and known higher order SM corrections as computed in 



ref. |31|, as well as the gluonic two-loop MSSM ones ||3j. (In the CMSSM two-loop gluino 
corrections are typically subdominant because the colored superpartners tend to be rather 
heavy.) These updates and improvements lead, however, to fairly minor changes in the 
overall distribution of most probable CMSSM parameter regions relative to refs. [l^, in 
agreement with ref. ||T2[| . On the other hand, as we have mentioned, the recent downwards 
shift in the SM value of BR[B — > Xs'^) has caused a corresponding big change in the pdf 
distribution in the CMSSM parameters |1£]. 



Finally, points that do not fulfil the conditions of radiative EWSB and/or give non- 
physical (tachyonic) solutions are discarded. We adopt the same convergence and mixing 



criteria as described in appendix A2 of ref. ||T^, while our sampling procedure is described 



in appendix Al of ref. [111] ]. We have the total of A'^ = 10 MC chains, with a merged number 
of samples 3 x 10^, and an acceptance rate of about 2%. More details of our numerical 



MCMC scan can be found in |11]. 



^We note that, in contrast to refs. |l^, |l^, in translating LEP-II bounds in the low jtia regime into 
our likelihood function we do not assume a priori that the CMSSM light Higgs scalar is SM-like. 
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Figure 1: The 1-dim relative probability densities for Xt (left panel) and Xb (right panel). Here 
and in all subsequent figures all parameters which are not shown have been marginalized (i.e., 
integrated) over. 



Roszkowski, Ruiz & Trotta (2006) 




Relative probability density 

0.2 0.4 0.6 0.8 1 



Figure 2: The 2-dim relative probability density in the plane of Xt and Xt- The probability color 
code shown here applies also to all subsequent figures showing 2-dim relative probability density. 
The inner (outer) solid contours delimit the regions of 68% and 95% total probability, respectively. 

3. Properties of the lightest Higgs boson in the CMSSM 

In this work we are particularly interested in the properties of the lightest neutral Higgs 
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boson. At the tree level, the Higgs sector of the MSSM is determined by tan/3 and niA- 
The (by far dominant) one-loop radiative corrections are generated by diagrams involving 
the top quark and its scalar partners, and, at large tan/3, also the bottom quark and 
its scalar partners. Both are proportional to their respective Yukawa couplings. The 
radiative corrections have been computed using several different methods. Full one-loop 
expressions are known |33, 34, Leading two-loop corrections have been computed 

using renormalization group [36| and two-loop effective potential methods 38 1, and in 
the Feynman-diagrammatic approach |3^]. Furthermore the tadpole corrections, needed 
to minimise the effective scalar potential, have been calculated at one loop 40 1 and the 
leading ones at two loops E3] . The remaining theoretical uncertainty in the light Higgs 



mass mfi has conservatively been estimated at < 3 GeV |43, 44 1. 

In computing the Higgs (and SUSY) mass spectrum we employ the code SOFT- 



SUSY v2.08 1 45 1, which implements radiative corrections in the modified Dimensional 



Reduction scheme, DR, based on the results of |35, 38, 41 1. For comparison, in Feyn- 
Higgs 1 46 1 the On-Shell scheme approach is adopted. Both are in good agreement [43|, i.e., 
within the above theoretical uncertainty. 

In the literature one often considers the cases of "maximal mixing" and "no mixing" 
(or "minimal mixing" ) to describe the impact on the Higgs sector of the off-diagonal terms 
rUqXq (q = t, b) in the stop and sbottom mass matrices relative to their diagonal entries. 
In terms of 

Xt = At — ^ cot (3, Xh = Ab — ^tan f], (3-1) 

where At^b are the stop/sbottom trilinear soft parameters, the "no mixing" case, in particu- 
lar, corresponds to = and Xb = 0. In fig. || we show in the left panel the 1-dim relative 
probability densities for Xt and in the right panel those for Xb in the CMSSM. Their 2-dim 
relative probability density (marginalized over all other parameters) is shown in fig. |2|. One 
can see that both variables are typically negative, with \Xb\ ^ Very large negative 

values of Xb are predominantly caused by the fact that the relative probability density of 
tan/3 is strongly peaked at large values of ~ 52 11, |l^; compare fig. 2 of ||ll|. Such 
large values of both Xt and Xb do not, however, necessarily imply large mixings in the stop 
and sbottom sectors since in the corresponding mass matrices they are multiplied by their 
respective quark masses. In fact, in the CMSSM in the sbottom sector we find a nearly 
strict no-mixing limit while in the stop sector we find a spread of values (mild mixing) 
but again with the peak in the probability distribution close to no mixing. Note that the 
ranges of Xt and Xb in figs. |^ and |2| are very different from the values proposed for general 
MSSM Higgs searches at the Tevatron [ 



The values of Xt and Xb determine to some extent the upper bound on rrih 



_2 



In ^ + 1 



mf 



■0 



12m 2 

t , 



+ (t^b,t 



(3.2) 



where g2 is the SU{2)l gauge coupling and m~ = {m~ + m~ )/2 is an average stop mass- 
squared, and analogously for the sbottoms. From fig. ^one can easily see that the bottom- 



sbottom contribution to (3.2) can be comparable to the top-stop one. 
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Roszkowski, Ruiz & Trotta (2006) 
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110 115 120 125 
m, (GeV) 



Figure 3: The 1-dim relative probability density for the light Higgs boson mass m/j. Magenta 
(dark shaded) and cyan (shaded) bands delimit 68% and 95% posterior probability (2 tails) ranges 
of rrih, respectively. 



Roszkowski, Ruiz & Trotta (2006) Roszkowski, Ruiz & Trotta (2006) 




Figure 4: The 2-dim relative probability densities in the plane of (mo,m/i) (left panel) and 
('^1/2 1 'Ti/i) (right panel). The inner (outer) solid contours delimit the regions of 68% and 95% total 
probability, respectively. 



In fig. ^ we display the 1-dim relative probability density of the light Higgs scalar 



-10- 



Roszkowski, Ruiz & Trotta (2006) Roszkowski, Ruiz & Trolta (2006) Roazkowski, Ruiz Sc TroUa (S006) 




Figure 5: The 1-dim relative probability densities for the masses of H, A and H^, respectively. 
Magenta (dark shaded) and cyan (shaded) bands delimit mass ranges of 68% and 95% of posterior 
probability (2 tails), respectively. 



mass. It is clearly well confined, with the ranges of posterior probability given by 

115.4 GeV <mh< 120.4 GeV (68% region), 

(3.3) 

112.5 GeV <mh< 121.9 GeV (95% region). 

The finite tail on the l.h.s. of the relative probability density in fig. below the 
final LEP-II lower bound of 114.4 GeV (95% CL) is a consequence of the fact that our 
likelihood function does not simply cut off points with m/j below some arbitrary CL, but 
instead it assigns to them a lower probability, as described above. On the other hand, 
the sharp drop-off on the r.h.s. of the relative probability density is mostly caused by 
the assumed upper bound on niQ < 4 TeV. This is shown in the left panel of fig. The 
upper bound on nih increases with nij and (3.2), whose largest values in turn depend 



on the maximum allowed value of mo, as we could already see in fig. 5 of |11]. With the 
new SM value for BR{B — > Xs'^) in the current analysis, the dependence on the prior 
range of mo has, unfortunately, become even stronger. For instance, adopting a much 
more generous upper limit mo < 8 TeV would lead to changing the ranges |3.3| to roughly 
120.4 GeV <mh< 124.4 GeV (68% CL) and 115.4 GeV < m^ < 125.6 GeV (95% CL). 
We will come back to this issue when we discuss light Higgs detection prospects at the 
Tevatron. 

On the other hand, the upper bound on m^ does not depend on the upper limit 
imposed in the prior for mj^/2) can be seen in the right panel of fig. ^ (compare also 
fig. 2 in ref. |0] ) . Basically, for very large tan /3 > 60 the bottom quark Yukawa running 
coupling becomes non-perturbative below the unification scale and it is no longer possible 
to find consistent mass spectra using the RGEs. This upper bound on tan/3 limits from 
above the values of mi/2 that can still be consistent with ^^^"^ l^-*-!- 



Fig. ^ confirms the well-known fact that in the CMSSM the largest values of rrih are 
typically much lower than in the general MSSM where stop and sbottom masses can to a 
large extent be treated as free parameters. The shape of the relative probability density 
also agrees rather well with ref. [12|. 
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Figure 6: The 2-dim relative probability density in the plane of (to^, m^) (left panel), (to/j, tan /3) 
(middle panel) and (to^i, tan/3) (right panel). 

The other Higgs bosons are typically somewhat, but not necessarily much, heavier. 
This can be seen in fig. ^, where we show the 1-dim relative probability densities for the 
masses of H, A and H^, respectively. Note that the shapes of their relative probabil- 
ity densities are almost identical since the masses of the three Higgs bosons are nearly 
degenerate. Their posterior probability regions are given by 

0.4 TeV < mH,A,H± < 2-5 TeV (68% region), 

0.2 TeV < mH,A,H± < 3-6 TeV (95% region). ^^'^'^ 

It is interesting to see a 2-dim relative probability density in the (rrih, plane. This 
is presented in the left panel of fig. ^. One can see that the pseudoscalar mass tua always 
remains somewhat larger than ruh, but that most of its probability density is actually 
concentrated at fairly low values. (Compare also fig. |5|.) For rrih < 115 GeV (118 GeV), 
at 68% CL we find ttia < 0.8 TeV (1.5 TeV). Experimental limits on Higgs boson searches 
in the MSSM are often presented in the plane spanned by the mass of the Higgs boson 
and tan /?. Our CMSSM results for the 2-dim relative probability density are presented 
in fig. ^ in the plane of (m/j,tan/?) (middle panel) and (m^,tan/3) (right panel). It is 
clear from fig. |6| that in the CMSSM uia ^ mz (the decoupling regime) but we find tua 
predominantly in the regime of a few hundred GeV, which is on the borderline of a "mild 
decoupling" regime. This will affect some relevant Higgs couplings, as we will see shortly. 

At this point we want to digress to emphasize the difference between posterior proba- 
bility (Bayesian statistics) and the mean quality-of-fit statistics. In the absence of strong 
constraints from data, the two can produce quite different distributions, which ought to 
be interpreted carefully since their meaning is different. (See ref. |11] for more details and 



a thorough discussion.) We illustrate this in fig. ^, for the case of Higgs masses where 
we show the distribution of the mean quality-of-fit (dashed line) in addition to the rela- 
tive probability density. Next, in fig. |^ we present the mean quality-of-fit in the planes of 
{rrih, tua) (left panel), {ruh, tan 13) (middle panel) and (m^, tan /?) (right panel). This figure 
should be compared with fig. |^. It is clear that the largest values of the mean quality-of-fit 
show preference for a smaller ruh, with the peak of the distribution at 114 GeV (below 



- 12 - 



Roszkowski, Ruiz & Trotla (2006) 



Roazkowskl, Ruiz & Trotli 




110 115 120 125 

m, (GeV) 




Roszkowski. Ruiz & Trotta (2006) 



Roszkowski, Ruiz & Tralta (SOOS) 



0.8 



0.6 



D. 0.4 



K 0.2 




1 2 3 
(TeV) 




1 2 3 
m„. (TeV) 



Figure 7: The same as in figs. || and || but with the mean quality-of-fit (dashed line) added for 
comparison. 
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Figure 8: The mean quahty-of-fit (dashed line) in the planes of (m/i, toa) (left panel), (m^, tan/3) 
(middle panel) and (m^,tan/3) (right panel). The 68% and 95% posterior probability contours 
from fig. [g] are also shown for comparison. 



the LEP-II bound of 114.4 GeV). The mean quality-of-fit also favors significantly smaller 
ruA < 1 TeV (as well as the other heavy Higgs masses). We note that, according to the 
mean quality-of-fit statistics, the favored region of parameter space lies at smaller masses 
than the 68% range of posterior probability, as can be seen in fig. |^. Such a discrepancy 
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Figure 9: The 2-dim relative probability density in the plane of tan /3 and the decay width T{h 
&&)cMSSM (left panel) and r(/i — > t+t^ )cmssm (middle panel) in the CMSSM, both normahzed to 
the SM values. In the right panel the same quantity is plotted for the ratio of the two widths. The 
inner (outer) solid contours encompass regions of the 68% (95%) total probability, respectively. 



between the two statistical measures can only be resolved with better data. We also note 
that we have found a handful of points (about 100) with C| <C 1 exhibiting a very good 
quality-of-fit at very small values of nih ~ 90 GeV. Since their statistical weight is insignif- 
icant (compared to some 3 x 10^ samples in our chain) we do not display their quality-of-fit 
in fig. g 
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Figure 10: The same as in fig. ^ but with tan/3 replaced by ruA- 
Next, we discuss some relevant couplings. We compute them with the help of Feyn- 



Higgs v2.5.1 |46] via the corresponding decay widths. For this purpose we translate all 
the Higgs and superpartner mass spectra and other relevant parameters from the output 
of SOFTSUSY. The decay widths of the neutral Higgs bosons are evaluated with the full 
rotation to on-shell Higgs bosons, i.e., beyond the effective coupling approximation. That 
means that the mixing matrix is constructed through the Z-factors resulting from the 
Higgs boson wave function normalization In the case of the vector bosons V = Z,W 
the coupling ratio is g{hVV)MSSM/ g{hVV)sM = sin(/3 — Ogfj), where a^^ denotes the ef- 
fective (i.e., radiatively-corrected) mixing angle in the Higgs scalar sector. The ratio is 
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typically strongly suppressed in the general MSSM. In contrast, in the CMSSM it is close 
to 1 as we have already mentioned when discussing C,f^ in the previous section.^ This is the 
case in both the decoupling and the mild decoupling regimes. 

The light Higgs couplings-squared to the third generation down-type fermions show 
a more complex behavior, as shown in figs. ^ and In the (C)MSSM'^ the light Higgs 
coupling to bottoms 5(/i6&)cmssm and to taus g(/iT"''r~)cMSSM, normalized to their SM 
value, g{hff)sM = g^nb^r/'^f^w-: are given by — sinaeg/ cos/3 = sin(/3 — aeff) — tan/3cos(/3 — 
aefr). At tree level both ratios are equal to 1. We note two effects here. First, for not too 
large tan /3 < 45 (for which in the CMSSM there is strong preference for > 1 TeV) we 
find that the ratio TcMSSM/rsM becomes close to 1.01. This is probably again just a result 
of an approximation used in FeynHiggs, as in the case of the VV mode. At larger tan (3 
the second term starts playing a bigger role. As ruA decreases to below some 0.8 TeV, 
corresponding to large tan (3 > 50, both coupling-squared ratios grow rather fast. The 
enhancements in this region can be seen in the left and middle panels of figs. |9| and |l^. At 
the 2(7 level, we find that the widths are increased relative to their SM counterparts by a 
factor of up to some 2.5. 

The second effect on the couplings is caused by radiative corrections from sbottom- 
gluino and stop-higgsino loops to the tree-level relation between the bottom mass and 
its Yukawa coupling |47]. At large tan/3 this leads to modifying the coupling g{hbb) p9|, 



while the analogous coupling to taus is not affected. (Implications of this effect for Tevatron 
Higgs searches have recently been discussed in ref. [^.) As a result, in the CMSSM at 
large tan /3 > 50, in a sizable number of cases the quantity g'^{hbb) cmssm, while remaining 
dominant, will show a small decrease relative to 5'^(/it+t~)cmssm- This feature, which is 



displayed in the right panels of figs, u and 10, will give one some chance of producing a 



somewhat increased number of taus in light Higgs decays at the Tevatron, as we will see 
shortly. 

We will now examine the impact of the above properties of the light Higgs couplings 
on its decays. The varying of the couplings to bottoms and staus is reflected in the total 
and partial decay widths of h to bb and r^r^ for which relative probability densities are 



presented in fig. 11. For comparison, we show the same quantities for the SM Higgs with 
the same mass. Somewhat larger widths in the bb and t~^t~ modes are a result of both 
decay channels being enhanced at large tan f5. On the other hand, we have checked that 
the decay width for h — > W*W (followed by vector boson decays into light fermions) shows 
no deviation from the SM. 

The ensuing effect on the branching ratios is shown in figs. |l^-0. We remind the 
reader that, in the range of mass predicted in the CMSSM, the SM-like light Higgs boson 
decays predominantly into bb pairs (~ 90%), followed by r+r" pairs (~ 9%), although at 
small tan /3 the h WW* branching ratio grows quickly with ruh and at ruh ~ 120 GeV 



^We do find a tiny increase, at the level of 1.5%, which is probably an artefact of the approximations 
used in FeynHiggs for computing a^s- 

'^From now on we will denote several MSSM variables with the subscript "CMSSM" , in order to emphasize 
the fact that their numerical values presented here are specific to the CMSSM, rather than to the general 
MSSM. 
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Figure 11: The 1-dim relative probability densities for the light Higgs total decay width (left 
panel) and partial decay widths to bb (middle panel) and t~^t~ (right panel). For comparison, the 
widths of the SM Higgs boson with the same mass are also shown. 




Figure 12: The 1-dim relative probability density of the branching ratios BR{h bb) (left 
panel), BR(h t^t^) (middle panel) and BR{h WW*) (right panel), all normalized to their 
SM counterparts. 
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Figure 13: The 2-dim relative probability density of the branching ratios BR{h bb) (left 
panel), BR{h — > t+t~) (middle panel) and BR{h WW*) (right panel), all normalized to their 
SM counterparts, as a function of m/j. 
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Figure 14: The 2-dini relative probability density of the branching ratios BR{h — > bb) (left 
panel), BR{h t+t~) (middle panel) and BR{h WW*) (right panel), all normalized to their 
SM counterparts, as a function of tan/3. 



it can exceed some 10% (at the expense of the above two channels). The dominance of the 
bb mode, however, is so large that, despite some decrease of the coupling hbb relative to the 
one of hT~^T~ at large tan /? in some parts of the CMSSM parameter space, the branching 
ratio into bb remains basically unaffected. On the other hand it is the subdominant modes 
r+r~ and WW* that, at large tan/3, can experience either a relative increase or decrease, 
respectively. These effects can be seen in fig. where we display 1-dim relative pdf's of 
the branching ratios for h ^ bb (left panel), h — > r^r^ (middle panel) and h — > WW* 
(right panel), all normalized to the analogous quantities for the SM Higgs boson. The 
corresponding 2-dim relative probability densities are shown in figs. 13 and|l^ in the planes 
spanned by the above SM-normalized branching ratios and nih and tan/?, respectively. 
Note a small increase in the number of produced r-leptons (which is rather small to start 
with), which may help in Higgs searches in that important decay channel. 



In conclusion, in the CMSSM with fiat prior ranges as given at the beginning of sec. 2.2 



(most importantly mo < 4 TeV) and with observables as in tables |l] and the mass of 
the light Higgs lies predominantly in the range shown in (3^) and fig. |^, while the other 
Higgs bosons are typically somewhat, but not necessarily much, heavier. The light Higgs 
coupling to VV remains basically SM-like, while the ones to bb and r+r" do show some 
increase relative to the SM ones at large tan /?. This has some effect on on light Higgs 
decays and, in the case of the bottoms, also on its production, as we shall see in the next 
section. 



4. Light Higgs production and decay 

We will now assess the discovery prospects of the light CMSSM Higgs boson at the Tevatron. 
To this end we will consider the following production and decay processes: 

• vector boson bremsstrahlung: V* Vh (where V = W, Z), followed by /i — > bb, 

h T^T^ or, in the case of Wh, also h WW*; 

• gluon— gluon fusion: gg — > h, followed by either h — > WW* or /i — > r+r^; 
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• associated bottom production: hb{b), with a b quark tagged on a hard spectrum 
(pt > 15 GeV and r] < 2.5), followed by either /i ^ 66 or /i — > t^t^; 

• inclusive production: pp h, followed hy h ^ t^t^ . 

We have also considered VV fusion and tih Higgs production processes, but in the 
CMSSM these are subdominant. 

Some additional comments about the processes that we consider are in order. The 
vector boson bremsstrahlung process is determined by the effective coupling g[hVV) ~ 
sin(/3 — Oeff). On the other hand, the other three processes are to a large extent determined 
by the behavior of the effective coupling g{hbb) which, as we have seen above, for tan /? > 50 
and mA % 1 TeV, can markedly deviate from the SM value. 

In the gluon-gluon fusion process gg ^ h we include diagrams with top or bottom 
quark (and their superpartner) lines in the loop. In the associated bottom production 
process hb{b) we compute the cross section of bg bh. An alternative, and effectively 
equivalent, way would be to consider the process gg — > hbb with the momentum of one of 
the bottoms integrated out [^]. The inclusive process pp ^ h can likewise be computed 
in two ways. At NLO in the four-flavor scheme one can add the (dominant) process 
gg hbb and the (subdominant) one qq hbb, and integrate out the momenta of the 
bottoms. Alternatively, one obtains very similar results by computing the process 66 — > /i 
at NNLO in the five-flavor scheme [^]. Here we follow the latter approach. Relative 
strengths of the above processes depend on several parameters, especially on tan /? (when 
large) but typically, in the MSSM in the regime of light Higgs mass and large tan /3, the 
gluon-gluon fusion process is dominant and associated bottom production is a factor of 
a few smaller but otherwise comparable. For a comprehensive review of Higgs properties 
and collider search prospects, see, for instance, ref. Q. 

Since, as we have seen above, the light Higgs boson is SM-like, it will be convenient 
to normalize our results to the corresponding processes involving the SM Higgs boson with 
the same mass. We expect most ratios to be close to 1 and in fact it is some possible 
departures from the SM case that we will attempt to identify. 

We compute the light Higgs production cross sections in the CMSSM, normalized 
to their SM counterparts, (Tcmssm/csm, with the help of FeynHiggs v2.5.1 [46|. The 



package implements the Higgs production cross sections at the Tevatron and the LHC, both 
evaluated in the effective coupling approximation using the SM cross sections provided in 



ref. [^] . The calculation of the branching ratios is based on ref. 1 52 ] . The code also includes 
SUSY corrections to Higgs couplings to the bottom quarks, which can be substantial at 
large tan /?, as discussed earlier. 

As above, we will follow the procedure developed in ||ll[ in presenting our results in 
terms of relative posterior pdf, here simply called relative probability density, for various 



variables. First, in fig. |l5|we display relative probability densities for <tcmssm/csm for the 
processes of primary interest at the Tevatron. All the ratios are close to 1 but only in 
the case of V* — > Vh {V = Z, W) is a pdf very strongly peaked very close to 1. The pdf 
for the gluon-gluon fusion SM-normalized cross section is peaked around 1.1, with rather 
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Figure 15: The l-dim relative probability density for light Higgs production cross sections at the 
Tevatron (normalized to the SM case) for V* Vh, where V = Z,W (upper left panel), gg ^ h 
(upper right panel), hb{b) (lower left panel) and pp h (lower right panel). 

little variation. This is a reflection of the behavior of the (radiatively corrected) coupling 
g{hbb)cMSSM] compare the left panels of figs. |^ and For the remaining two processes we 
observe some more variation, and increase relative to the SM, than in gg — > h. Actually, 
because of the way these processes are computed (as described above), their pdf's will in 
most cases be very similar. We nevertheless present them both for completeness. 

In order to display the behavior of the production cross sections in more detail, we 



present in fig. |16| a 2-dim relative probability density of the SM-normalized cross sections 
as a function of rrih, in fig. ^ of tan/3 and in fig. |l^ of the pseudoscalar Higgs mass niA- 
As regards the gluon-gluon fusion process, the bottom quark exchange contribution to the 
cross section is subdominant relative to the top quark one (by a factor of a few). This 
explains why there is much less variation in the corresponding pdf than for hb{b) and 
pp — > h. On the other hand, the enhancement of the coupling g{hbb)cMSSM at tan/3 > 50 
and niA 1 TeV, cause a slight increase of a few per cent in the cross section relative to 
the SM. Otherwise, unsurprisingly, the pdf's for the three processes mirror the behavior of 
the g{hbb)cMSSM coupling and we include them for completeness. To finish our discussion 



of light Higgs production, we present in fig. |19| the ranges in the CMSSM of the SM 



normalized cross sections for the above two processes in the plane of tan/3 and m/j, while 



in fig. 20 the same is shown in the plane of tan /? and uia- 



We now combine the above results for the light Higgs production cross sections and 



decay branching ratios at the Tevatron. In fig. 21 we show l-dim relative probability 
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Figure 16: The 2-dini relative probability density of light Higgs production cross sections at the 
Tevatron (normalized to the SM case) as a function of its mass m/j. 

densities for SM-normalized light Higgs production cross section times decay branching 
ratio {a x BR)cmssm/{o' x BR)sm, while the ratio's dependence on ruh is displayed in 



fig. |2^, on tan (3 in fig. ^3], and on itla in fig. In fig. 25 we show a distribution of values 
of the above product in the plane spanned by rrih and tan /?, while in fig. ^ the same 
quantities are shown in the usual plane of and tan/3. As before, all parameters other 
than the ones shown in each figure have been marginalized over. 

The emerging picture is rather clear. As expected, in the CMSSM, for all the considered 
processes, we generally find very similar light Higgs search prospects as for the SM Higgs 
boson with the same mass. We note, however, some differences, which may help optimize 
search strategies. To start with, in the vector boson bremsstrahlung process V* — > Vh 
{V = Z, W), h — > bb, T'^T' modes are almost indistinguishable from the SM Higgs case with 
the same mass. This is caused by the fact that in the CMSSM the coupling g{hVV)cMSSM 
is very close to its SM value. We note, however, that we do find a slight enhancement 
of the bb and t~^t~ final states, which may be of some help in these important search 
channels. The same is of course true for all the other processes considered here. On 
the other hand, we have found in the CMSSM parameter space some deviations in the 
g{hbb)cMSSM coupling from the SM value, which may significantly change production cross 
section of all the modes except V* Vh. On the other hand, the h — > WW* mode is 
typically reduced by up to some 5% within the 68% posterior probability region. Within 
the 95% posterior probability region, it can actually be significantly reduced, especially in 



-20- 



Roszkowski, Ruiz S Trotta (2006) Roszkowski. Ruiz & Trotta (2006) 




tanp tanp 



Figure 17: The 2-dini relative probability density for light Higgs production cross sections at the 
Tcvatron (normalized to the SM case) as a function of tan/?. 



the region of tan /5 > 50 and uia ^ 1 TeV. 

We emphasize that the results presented in figs. have been derived in the frame- 



work of the CMSSM. This should be kept in mind when comparing them with existing 
experimental Higgs search limits. Many of them have been set for specific (e.g., maximal 
and no mixing) scenarios and/or choices of parameters in the general MSSM, which are ba- 
sically never realized in the CMSSM. Conversely, it would be helpful to add to experimental 
results limits applicable to the CMSSM. 



5. Summary and conclusions 

In this work we have performed a global scan of the CMSSM parameter space by employing 
a powerful MCMC technique. We have then analyzed our results for light Higgs properties 
and discovery prospects at the Tevatron, mostly in terms of Bayesian statistics, although 
we have demonstrated that an alternative mean quality-of-fit analysis can lead to rather 
different results. In particular, while the former favors the light Higgs mass range above 
the final LEP-II 95% CL, the latter points more towards values below it for a large part. 

The couplings of the light Higgs of the CMSSM to vector bosons and bottoms and 
taus are basically very similar to those of the SM Higgs boson with the same mass. Small 
enhancements, at the level of a few per cent, have been found in most of the CMSSM 
parameter space, although at large tan/? > 50 and niA < 1 TeV (the preferred region). 
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Figure 18: The 2-dim relative probability density for light Higgs production cross sections at the 
Tevatron (normalized to the SM case) as a function of tua- 
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Figure 19: Values in the (m/i, tan f3) plane of light Higgs production cross sections at the Tevatron 
(normalized to the SM case) for the processes hb{b) and pp h. For the processes V* — > Vh, 
where V — Z, W, and gg —> h the values are almost independent of tan/3; compare fig. 



the differences can be substantial. Our intention was to provide experimentalists involved 
in Higgs searches at the Tevatron detailed information about detection prospects of the 
CMSSM light Higgs boson. Despite the fact that the pdf 's for associate bottom production 
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Figure 20: Values in tlie {niA , tan /3) plane of light Higgs production cross sections at the Tevatron 
(normalized to the SM case) for the processes hb{b) and pp — > h. For the processes V* Vh, where 
V = Z, W , and gg h the values are almost independent of tan/? and m^; compare figs. |l7| and [l^. 



and inclusive Higgs production modes are basically indistinguishable, we have displayed 
them separately for the sake of completeness and convenience. 

At the Tevatron, the sensitivity to a SM-like Higgs boson in the mass range of up to 
some 122 GeV (compare ( |3.3| )) seems excellent. According to ref. |^3[, with about 2fb^^ of 
integrated luminosity per experiment (expected by the end of 2006), a 95% CL exclusion 
limit can be set for the whole 95% posterior probability light Higgs mass range given *** 
in ( ^.3| ) ****. We stress again that this conclusion depends on our assumed ranges of 
flat priors, especially on niQ < 4 TeV, as discussed in the text. (Extending the range 
of niQ up to 8 TeV, and accordingly nih up to < 125.6 GeV (95% CL), would require 
instead about 2.5 fb~^ of integrated luminosity per experiment, which is again well within 
Tevatron's reach.) While keeping this in mind, we still find it remarkable that negative 
Higgs searches at the Tevatron should allow one to make definitive conclusions about the 
ranges of CMSSM parameters, in particular mo, which extend well beyond the reach of 
even the LHC in direct searches for superpartners. 

Should a signal (hopefully) start being seen in this mass range, a 3a evidence {5a 
discovery) can be claimed with about 4fb~^ (12 fb^^) per experiment. On the other hand, 
with about 8fb~^ of integrated luminosity ultimately expected per experiment, a 5c7 discov- 
ery will be possible, should the light Higgs mass be around 115 GeV. While such low mass 
is just below the 68% region of posterior probability according to Bayesian statistics, it is 
actually favored by an alternative mean quality-of-fit analysis. In other words, according 
to this measure, light Higgs search appears more promising than in the more conservative 
Bayesian probability scenario. A conclusive search for the CMSSM light Higgs boson at 
the Tevatron seems therefore fully feasible. 
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Figure 21: The 1-dim relative probability (iensity for light Higgs production cross sections times 
decay branching ratios at the Tevatron (normalized to the SM case). 
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Figure 22: The 2-dim relative probability density of light Higgs production cross sections times 
decay branching ratios at the Tevatron (normalized to the SM case) as a function of its mass m/j. 
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